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Using the Mathisson-Papapetrou-Dixon (MPD) equations, we investigate the trajectories of a 
spinning particle starting near r^^' in a Kerr field and moving with the velocity close to the velocity 
of light (r^^-* is the Boyer-Lindquist radial coordinate of the counter-rotation circular photon orbits). 

First, as a partial case of these trajectories, we consider the equatorial circular orbit with r = r^^'. 
This orbit is described by the solution that is common for the rigorous MPD equations and their 
linear spin approximation. Then different cases of the nonequatorial motions are computed and 
illustrated by the typical figures. All these orbits exhibit the effects of the significant gravitational 
repulsion that are caused by the spin-gravity interaction. Possible applications in astrophysics are 
discussed. 

PACS numbers: 04.20.-q, 95.30.Sf 



I. INTRODUCTION 

The geodesies in a Kerr metric are considered in the 
classical books on general relativity [1-3]. Some recent 
papers are devoted to more detailed study of geodesies 
on Kerr's black hole with the aim to elucidate the mech- 
anism of jet formation [4], and to analyze the possibil- 
ity of particle acceleration to arbitrary high energy [5]. 
The complete sets of analytic solutions of the geodesic 
equation in axially symmetric space-time are given in [6] . 
However, the description of particle motion by geodesies 
is restricted to a spinless particle. The motion of a 
spinning test particle is described by the Mathisson- 
Papapetrou-Dixon equations [7-9] : 



D 

ds 



ds 



ds J 



ds 



ds 



0, 



(1) 
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where = dx^ /ds is the particle's 4- velocity, S'^'^ is 
the tensor of spin, m and D/ds are, respectively, the 
mass and the covariant derivative with respect to the 
particle's proper time s; R^p^, is the Riemann curvature 
tensor (units c—G—1 are used) . It is necessary to add 
a supplementary condition to Eqs. (1), (2) in order to 
choose an appropriate trajectory of the particle's center 
of mass. Most often the conditions 



and 



are used, where 



S-^^u. = 



S^^'Py = 



= mu -\- ux 



DS' 



ds 
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(4) 



(5) 



is the 4-momentum. In practice, the condition for a spin- 
ning test particle 



£ < 1 



(6) 



must be taken into account [10], where jS'oj = const is 
the absolute value of spin, r is the coordinate distance of 
the particle from the massive body. 

After [7-9], Eqs. (1), (2) were obtained in many pa- 
pers by different approaches [11]. Also, this subject is of 
importance in some recent publications [12]. 

In general, the solutions of Eqs. (1), (2) under con- 
ditions (3) and (4) are different. However, in the post- 
Newtonian approximation these solutions coincide with 
high accuracy [13], just as in the case if the spin ef- 
fects can be described by a convergent in spin series, as 
some corrections to the corresponding geodesic expres- 
sions [14]. Therefore, instead of rigorous MPD Eqs. (1) 
their linear spin approximation 



m—u 

ds 



(7) 



is often considered. In this approximation condition (4) 
coincides with (3). 

The effect of spin on the particle's motion in Kerr's 
field has been studied since the 1970s [10,15,16]. In the 
past 10-12 years this subject gives rise to renewed inter- 
est [17-22], particularly in the context of investigations 
of the possible chaotic motions [17,19]. Also, these refer- 
ences provide a good introduction concerning the MPD 
equations. 

The purpose of this paper is to investigate more care- 
fully the world lines and trajectories of a spinning particle 
moving relative to a Kerr source with the velocity close 
to the velocity of light. We focus on the circular and 
close to circular orbits, because just on these orbits one 
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can expect the significant effects of tlie spin-gravity in- 
teraction [10,23-25]. Indeed, tliese orbits are of interest 
in the context of investigations of the nongeodesic syn- 
chrotron electromagnetic radiation of highly relativistic 
protons and electrons near black holes. Besides, it is 
known that the highly relativistic circular orbits of a spin- 
less particle are of importance in the classification of all 
possible geodesic orbits in a Kerr spacetime. Naturally, 
the circular highly relativistic orbits of a spinning parti- 
cle are of importance in the classification of all possible 
significantly nongeodesic orbits in this spacetime as well. 
Also, these orbits are exclusive in the sense that they are 
described by the strict analytical solutions of the MPD 
equations. The main features of the spin-gravity interac- 
tion that are revealed on the circular and close to circular 
orbits will be a good reference to further investigations 
of most general motions of a spinning particle in a Kerr 
spacetime. 

We stress that MPD equations are the classical limit 
of the general relativistic Dirac equation [26], new results 
in this context are presented in [27]. Therefore, highly 
relativistic solutions of the MPD equations stimulate the 
corresponding investigations of the fermion's interaction 
with the strong gravitational field. 

The paper is organized as follows. Sec. 2 deals with 
the relationships following from the MPD equations for 
the highly relativistic equatorial circular orbits of a spin- 
ning particle in Kerr's field, in the Boyer-Lindquist coor- 
dinates. The linear spin MPD equations for any motions 
of a spinning particle are considered in Sec. 3. The re- 
sults of computer integration of these equations for some 
significantly nongeodesic motions are presented in Sec. 
4. We conclude in Sec. 5. 

Following [3], in this paper notes the radial co- 
ordinate of the photon circular orbits in the case of the 
counter-rotation. 



where a dot denotes differentiation with respect to the 
proper time s, and square brackets denote antisym- 
metrization of indices; F^^ are the Christoffel symbols. 
Eq. (7) in terms of 5, is 



-.{u^Ki, + 2uiR^,^)e'''Si = 0. (10) 



In many papers the 4-vector of spin s\ is considered, 
where 

The following relationship holds: 

Si = UiS4 - U4Si. 



Let us consider Eqs. (9), (10) for the Kerr metric 
using the Boyer-Lindquist coordinates = r, = 
6, = x^ = t. Then the nonzero components of 
Qixv are 
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(11) 



where 



II. HIGHLY RELATIVISTIC EQUATORIAL 
CIRCULAR ORBITS OF SPINNING PARTICLES 
IN A KERR FIELD ACCORDING TO 
APPROXIMATE AND RIGOROUS MPD 
EQUATIONS 

In practical calculations it is convenient to represent 
the MPD equations through the spin 3- vector Si, instead 
of the 4-tensor S^" , where by definition 



Si = ^V^SikiS''^ 



(8) 



where g is the determinant of , Siki is the Levi-Civita 
symbol (here and in the following, latin indices run 1, 2, 
3, and greek indices 1, 2, 3, 4, unless otherwise specified). 
Then Eqs. (2) can be written as [23] 



U4Si - iiiSi + 2(u[4Uj] - u"upT''^^^UiT^)Skvf' 



+ 2SnVl^^Ui^u^ = 0, 



(9) 



p2 = ^2 ^2 ^.Qg2 A = r^-2Mr-^a^ < 6i < tt. 

[In the following we shall put a > 0, without any loss in 
generality; the metric signature is +.] It is easy 

to check that three equations from (9) have a partial 
solution with 6 = 7r/2, Si = Sr = 0, S3 = S,p = and 
the relationship for the nonzero component of the spin 
3- vector S2 = Se = is 



S2 = ruiSo, 



(12) 



where So is the constant of integration. The physical 
meaning of this constant is the same as in the general 
integral of the MPD equations [15] 



Q2 _ 1 C C/ii' 
"^0 — 2 '"'^ ■ 



(13) 



We stress that relationship (12) is valid for any equato- 
rial motions (6 = 7r/2), with the spin orthogonal to the 
motion plane {SiU^ = 0). 
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The possible equatorial orbits of a spinning particle are 
described by Eq. (10). First, we shall consider the case 
of the circular orbits with 

=0, = 0, = const ^0, = const =^ 0. 

(14) 

Investigating the conditions of existence of the equato- 
rial circular orbits for a spinning particle in Kerr's field 
we use Eqs. (10), (12), and u^u^ = 1. It is known that 
from the geodesic equations in this field, the algebraic re- 
lationship that follows determines the dependence of the 
velocity of a spinless particle on the radial coordinate r 
of the equatorial circular orbit. Similarly, from the first 
equation of set (10) using Eq. (14) we obtain the rela- 
tionship for the equatorial circular orbits of a spinning 
particle in Kerr's field as follows: 



A^S2 



x[{u^YR\,^+u'u\R\,,~Ri,,)-{uyRl,,]=Q. (15) 

[By Eq. (14) other equations of set (10) are automatically 
satisfied.] Taking into account Eq. (12) and the explicit 
expressions for T^^^ and iJ^^g. (see, e.g., [18]) from Eq. 
(15) we obtain 



3,. 4 



.4\2 



x[a^{r^+a^){uY-{r^+2a^)u''u'' + a{uy]^0. (16) 

The 4-velocity component can be expressed through 
from the condition u^U/j, = 1 as follows: 



544 



+ ■ 



/ 534^^533544. 
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(17) 



(just the sign "+" at the radical in Eq. (17) ensures the 
positive value of u^). Inserting expression (17) into Eq. 
(16) and eliminating the radical by raising to the second 
power we get 

{u^f{r^[r{r - 3Mf - 4Ma^] + 12eMor^A - Qe^MV^A) 
+{u^f{-2Mr'^{r - 3M) + 6sMar^{r - 3M) 

- 9£^MV^(r - 2M)) + M2(r - Seaf = 0, (18) 

where, as in Eq. (6), e = \SQ\/mr. Without any loss in 
generality we put 5*0 > 0, then by Eq. (12) 5*2 > 0. 

So, the particle's angular velocity ip = on the cir- 
cular orbit with the radial coordinate r must satisfy Eq. 
(18). 

Let us show that Eq. (18) provides known solutions. 
In the partial case of a spinless particle {s = 0) from Eq. 
(18) we have 

3 2 ^ Mrjr - 3M) + 2MaVM^ 
^ ' r2[r(r-3M)2-4Ma2] ' ^ ' 



It follows from Eq. (19) that the velocity of such a par- 
ticle on the circular orbit is highly relativistic if the ex- 
pression r(r — 3M)2 — AMa? is close to 0. This fact is 
known from the analysis of the geodesic orbits in a Kerr 
field, as well as that just the roots of the equation 

r(r - 3M)2 - AMa^ = 

determine the values of r for the photon orbits. 

If e 7^ and the absolute value of the expression 
r^[r{r - ?,Mf - AMa?] in the factor of (u^)'* in Eq. 
(18) is much greater than 12eMar'^A - Qe^M^r'^A, then 
it is easy to verify that the corresponding roots of Eq. 
(18) describe the circular orbits with the angular velocity 
which is close to the angular velocity of the correspond- 
ing geodesic orbits due to e ^ 1. More exactly, in this 
case we have 



/Mr(r - 3M) + 2Ma-/M 



-(l+0(£)), (20) 



where the main term is equal to the square root of the 
right-hand side of Eq. (19). 

Now we point out a case of interest with e 7^ that is 
not described in the literature. Namely, it is not difficult 
to calculate that for e 7^ Eqs. (16), (18) have the solu- 
tions which describe the highly relativistic circular orbits 
with the values of r that is equal or close to 



„(-) 

ph 



I.e., 

to the radial coordinate of the counter-rotation photon 
circular orbits. For example, in the case of the maximum 
Kerr field (a — M) the orbits with r ~ 4A'/(1 + 5), where 
I (5 1 <C 1, are highly relativistic, both for positive and neg- 
ative 5. If 1^1 ^ e, according to Eqs. (16), (18) the values 
of on these orbits are determined by the expression 



1 



3M^6£(1 -h 4(5/3e) 



(l + 0(£)). 



(21) 



[The choice of the sign in Eq. (21), < 0, is dictated 
by the necessity to satisfy both Eq. (18) and (16); Eq. 
(18), as compared to (16), has additional roots because of 
the operation of raising to the second power.] Similarly 
as in the case of Eq. (20), it is easy to check that if 
e <C 1^1 it follows from Eqs. (16), (18) the expression 
for which in the main term coincides with the known 
analytic solution for the corresponding geodesic circular 
orbit. 

It follows from Eqs. (16), (18) at r = for any 
< a < M that 



u = —- 



3£0 



.(-) 

ph 



M) 



(l + 0(£)). 



(22) 



It is known from the geodesic equations that the values 

of Tp^'' increase monotonically from 3M at a = to 4M 
at a = M. 

Thus, according to Eq. (22) the expression for the 
angular velocity <^ = w'^ in the main term is proportional 
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to '^l\fs, whereas the angular velocity in Eq. (20) at 

r = r"~i^(\ + (5), < (5 <C 1 is proportional to Xj^fb. 
Further details appear below in this Sec. 

Using Eq. (22) wc can estimate the value of the 
Lorentz 7-factor, corresponding to the 4- velocity compo- 
nent , for different a. More exactly, we shall calculate 
the Lorentz 7-factor from the point of view of an observer 
which is at rest relative to a Kerr mass. According to 
the general expression for the 3-velocity components 
we write [1] 



It is known that the important characteristic of the 
particle's motion in the Kerr spacetime are its energy and 
angular momentum. Let us estimate the conserved values 
of the energy E and angular momentum J of a spinning 
particle on the above considered highly relativistic circu- 
lar orbits with Eqs. (21), (22). The expressions for these 
quantities are [16,28] 

DS^^ DS'^^ 1 

E = mui+gsiUfj, +giiU^ + -S'^" 9vi,^, (29) 



dt + ^dx' 

.944 



(23) 



and for the second power of the velocity absolute value 
bp we have 



(24) 



where 7ij is the 3-space metric tensor. The relationship 
between 7ij and g^v is as follows: 



lij = 9ij + 



.944 



(25) 



For the circular motions we have =0, v"^ = 0, and 
according to Eq. (23) 



^/944 V 



544 



(26) 

By Eqs. (24)-(26) and the condition m'^m^ = 1, for the 
7-factor we write 



1 



-533 + ^1+1. 
544 



(27) 



Inserting the value from Eq. (22) into Eq. (27) we 
find in the corresponding spin approximation 



7 




-1/4 



2M 



-1/2 



(28) 



It follows from Eq. (28) that 7^ >> 1, i.e., the under con- 
sideration circular motions are highly relativistic. Fig. 1 
shows the dependence 7/70 on /M, where 70 is the 
7-factor for Tp^^ = 3M, that is, at a = 0. 




FIG. 1. Ratio 7-factor at different a, from < o < M, to 
7-factor at a = vs. radial coordinate. 



J = -TOW3 - g^zu^—^^ dsAUfi—^ -j^S^" guz,^i- 

(30) 



7 Cf-J^^u, 7 

as as 
By Eqs. (21), (29) for r = 4M(1 + 6), \S\ < £ we obtain 

TO 



6£ 



(31) 



The energy of a spinless particle on the geodesic circular 
orbit with r = 4M(1 -|- (5), < (5 < 1 for a = M is equal 
to 



E. 



3 m 



geod 



(32) 



It follows from Eqs. (31), (32) that 



i^geodM'»l. (33) 

At the same time according to (31), (32) for < ^ <C £ 
we have 



Ei 



eod 



20s 



(34) 



That is, the values of energy of the spinning and spinless 
particles on the highly relativistic circular orbits with the 
same r — 4M(1 + S) in the maximal Kerr field can differ 
significantly. It is easy to show that similar situation 
takes place for all values < a < M with r = rp'^\l + 6). 
In addition, one can estimate for these circular orbits that 
according to Eq. (30) J%iJ J%od < 1- 

As a result, the relationships E^^^^J E'^^^^ <C 1 and 

J%^n|Jleod < 1 foUowing from Eqs. (29), (30) show 
clearly that the corresponding solutions of the MFD 
equations cannot be obtained in the framework of the 
analytic perturbation approach to the dynamics of a clas- 
sical spinning particle developed in [20-22]. 

We stress that in this Sec. above we have considered 
the new highly relativistic circular solutions of the ap- 
proximate MFD Eqs. (2), (7). Let us show that these 
solutions satisfy the rigorous MFD Eqs. (1), (2). Indeed, 
it follows from Eqs. (1) that their terms, which were ne- 
glected in Eqs. (7), in the case of the circular equatorial 
motions are presented in the first Eq. of set (1) only. In 
metric (11) these terms can be written as 

Ci(«3)4 + c2(u4)^ + C3(«3)V + C4(«3)2(«3)2 + 05^3(^*4)3, 

(35) 
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where 



ci = SoA^{3r^ + a^){r^-Ma% 



C2 = -SqA zr 

r' 



A'Pa / 2M\ 



r ) 



C3 = SQ^Xr'ir - 3M)(r^ - SMa^) + 4:M^a% 
r' 



2/9 = 



Si 



2/10 



_S2_ 

mlVP ' 



and 



a; = 



M' 



£0 = 



2/11 = 



mM 



mM2' 



(37) 



(38) 



(in contrast to e from (6), that depends on r, here eq is 
defined to be const). Then it follows from Eqs. (7), (9) 
the set of 11 first-order differential equations 

2/1 = 2/5, m = 2/6, m = 2/7, 2/4 = 2/8, 



^5 = ^1, ye = ^2, y7 = ^3, ^8 = ^4, 



C5 



-5o 



MA 



3Mr^ - AMa^ + 



(36) 



First, we point out that according to Eq. (36) in the 
case of a Schwarzschild's field for the circular orbit with 
r = 3M the all coefficients Cj are equal to 0. There- 
fore, in this case expression (35) is equal to identically, 
independently on the explicit expressions for u^, . It 
means that the highly relativistic circular orbit of a spin- 
ning particle with r = ZM in a Schwarzschild's field is 
a common strict solution both of the approximate MPD 
equations (2), (7) and the rigorous MPD equations (1), 
(2). Second, it is not difficult to check that applying Eqs. 
(17), (22) to expression (35) at a ^ yields in the main 
terms, i.e., within the accuracy of order e. 

Thus, in this Sec. we dealt with new partial solutions 
of the MPD equations in a Kerr spacetime with Eqs. 
(12), (14), (22). All highly relativistic orbits of a spin- 
ning particle described by these solutions are circular and 
located in the space region withr = rj^'(l + ^), 1^1 « 1. 
To study highly relativistic orbits beyond this region it is 
necessary to carry out the corresponding computer cal- 
culations. It is an aim of Sees. 3 and 4. 



III. EQUATIONS (7), (9) FOR ANY MOTIONS 
IN A KERR FIELD 

Now the point of interesting is the noncircular highly 
relativistic motions of a spinning particle that starts near 
Tp^' in a Kerr field. In particular, we shall consider the 
effect of the 3- vector Si inclination to the equatorial plane 
Q = tt/2 on the particle's trajectory. In this case Eqs. 
(9) cannot be integrated separately from Eqs. (7). For 
computer integration, it is necessary to write the explicit 
form of Eqs. (7), (9) in metric (11). It is convenient to 
use the dimensionless quantities j/j, where by definition 



2/1 



r 

M' 



2/2 = 



2/3 = <P, 2/4 = 



t 

M' 



m = ^5, yio = ^6, yii = ^7, 



(39) 



where a dot denotes differentiation with respect to x, and 
Aj {j = 1, 7) are some functions which depend on t/j. 
Because the expressions for Ai in general case of any a 
are too lengthy, here we write Ai for the much simpler 
case a = 0: 

^1 = ^« + {y^yl + 2/12/7 sin' y2 - 



+- 



2/1 sin 2/2 



(2/72/10 sin' 2/2-2/62/11), 



A2 = 2/7 cos 2/2 sin 2/2 



22/52/6 32/72/9 . 

3— sin2/2, 

2/1 2/1 



. 22/52/7 3y62/9 

^3 = 22/62/7 cot 2/2 + -3-^ , 

2/1 2/1 sm 2/2 



. 22/52/8 3y5 2 

A4 = —q 3 : q 

2/1 2/1 2/8 sm 2/2 



X (2/72/10 sin' y2 - y&Vii), 



. ^ , 2/62/10 , 2/72/11 

^5 = —ysygq -\ 1 

2/1 2/1 2/1 



-(2/5^9 + yeyio + 2/7yii)[^ig 



2/5 . o2/5 2 

— A4q - 3-2 g 
2/8 2/r 



2/5 



ye = Mu', 2/7 = Mu^ 



2 2-2 ,2/81,2/9^ 

-2/12/6 - 2/12/7 Sin 2/2 + — H A4, 

yt 2/8 
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^6 = -ymvQ 1 + —9 + 1 



yw 



+2/72/11 cot 2/2 H ^4 - (2/52/9 + 2/62/10 + 2/72/ii) 

2/8 



X [ylA2 - + 2t/52/6(2/i - q) 



2/8 



terms. Wc monitor errors of computing using the con- 
served quantities: the absolute value of spin, the energy 
and angular momentum (see Eqs. (13), (29), (30)). 

Figs. 2-9 correspond to the case of Schwarzschild's 
field. All plots start with the same initial values of the co- 
ordinates = r, = 9, = x^ = t, namely, 
at 3M, 90°, 0°, and correspondingly. We do not 
vary the initial values of the 4-velocity components 
and as well. More exactly, we put = and 



-2/i2/7C0S2/2sin?/2] 



1 



3M\/2 



'-1 + 



VT2 



(42) 



, 2/52/11 / 2 
Ar = — { —q+1 



2/1 \2/i 



-2/12/72/9 I ^ ~ ^ 1 ®™ ^2 - 2/62/11 cot 2/2 



-2/72/10 cos 2/2 sin 2/2 - (2/52/9 + 2/62/io + 2/72/ii) [^32/i sin^ 2/2 



* 2/1 2/7 . 2 , n / \ • 2 

-A4 sm 2/2 + 22/52/7(2/1 - 9) sm 2/2 

2/8 



+22/1 2/62/7 cos 2/2 sin 2/2] -h — A4 , (40) 

2/8 

where 



2/1 



g= 1 



According to Eqs. (13), (38) the expression for £q can be 
written as 

el=y^q{l + qyl) + ^q'{l + ylyl) 
2/8 2/1 2/8 



2 2 

, yilQ r^ , 2 2-2 \ 

+ 22-2 — (1 + 2/12/7 sm 2/2) 
2/12/8 Sin 2/2 



2g2 

+ —(2/52/62/92/10 + ynyrysyn + 2/62/72/io2/ii)- (41) 

% 

Eqs. (39)-(41) are aimed at computer integration. 



IV. NUMERICAL RESULTS 

Wc present here the results of computer integration of 
Eqs. (39) with (40), (41). All plots below are restricted to 
the domains of validity of the linear spin approximation. 
That is, in these domains the neglected terms of the rig- 
orous MPD equations are much less than the linear spin 



where expression (42) is the solution of Eq. (18) at a = 0, 
rigorous in eo • That is, the initial values of and u"^ are 
the same as for the equatorial circular orbit with r = 3M 
(it is easy to check that expression (42) coincides with 
(22) in the corresponding spin approximation). However, 
we vary the initial inclination angle of the spin 3-vector 
to the plane 6 = 90°, without change of the absolute 
value of spin (Figs. 2 6), and add the small perturbation 
by the radial velocity (Figs. 7-9). Without any loss in 
generality we put Si > 0. For comparison, we present the 
corresponding solutions of the geodesic equations that 
start with the same initial values of the coordinates and 
velocity as the solutions of the equations for a spinning 
particle. In all cases as a typical value we put £0 = 10~^. 

Figs. 2-9 let us compare the world lines and 
trajectories of the spinning and spinless particles in 
Schwarzschild's field. Figs. 2, 6, 7, 9 exhibit the sig- 
nificant repulsive effects of the spin-gravity interaction 
on the spinning particle. Due to the repulsive action 
the spinning particle falls on the horizon surface during 
longer time as compared to the spinless particle (Fig. 2). 
Moreover, according to Figs. 6, 9 the considerable space 
separation takes place within a short time, i.e., within 
the time of the spinless particle's fall on the horizon. 

The point of interest is the phenomenon when a spin- 
ning particle orbits below the equatorial plane for some 
revolutions. Figs. 3, 8 (we recall that according to the 
geodesic equations similar situation is impossible for a 
spinless particle). It is easy to calculate that for 5i < 
a spinning particle can orbit above the equatorial plane. 

We stress that all graphs in Figs. 2-9 are interrupted 
beyond the domain of validity of the linear spin approxi- 
mation. By Figs. 3, 4, within the time of this approxima- 
tion validity, the period of the ^^-oscillation coincides with 
the period of the particle's revolution by (p. Whereas on 
this interval the value of the spin component Si is const 
(Fig. 5), just as the components ^2 and S3 (the corre- 
sponding graphs are not presented here). We point out 
that this situation differs from the corresponding case of 
the circular motions of a spinning particle that are not 
highly relativistic. Indeed, then the nonzero radial spin 
component is not const but oscillates with the period of 
the particle's revolution by (p. Besides, in the last case 
the mean level of 9 coincides with 90°, whereas in Figs. 3, 
8 the mean values of 9 are above 90°. According to Figs. 
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3, 8 the amplitude of the 6'-oscillation increases with the 
inchnation angle. However, 6 — 90° is small even for the 
inclination angle that is equal to 90°. 

In the context of Figs. 3, 8 we point out an interesting 
analytical result following from Eqs. (39)-(41). Namely, 
it is not difficult to check that these equations are satis- 
fied if 

3 

yi = - ^, 2/2 = arccos(52, 2/3 = 0, j/4 = 0, 

1 - di 

2/5 = 0, 2/6 = 0, y, = ±-^{l^6rf/\l-Sl)-y\ 

S^Jb^l 



spin (the inclination angle is equal to 0°) in Kerr's field 
at a = M. All curves start from r = 4M with zero radial 
velocity and the tangential velocity which at eo = 10~^ is 
needed for the circular motion with r = 4M. This circu- 
lar motion is shown by the horizontal line, whereas other 
curves represent the noncircular motions at eo < 10^^ 
with the same particle's initial values of coordinates and 
velocity. According to Fig. 12 the circular orbit of a 
spinning particle with r = AM monotonically trend to 
the corresponding noncircular geodesic orbit if Eq trend 
to 0, i.e., the limiting transition stq ^ is correct. 



ys 



2 
2/1 



-1/2 



3(5, 



y9 = ±-9^(l-5i)-'/'(l-<5i) 



-1/2 



2/10 



35i,. 



5i)-i(l-52)i/2, y^^=o, (43) 



where Si and 82 are some small constant values such that 

< (5i < 1, < \S2\ < 1. According to Eq. (41) the 
relationship between 61, 62 and s from Eq. (6) is as 
follows: 



(44) 



By notation (37) we conclude that the partial solution of 
Eqs. (39)-(40), that is presented in Eq. (43), describes 
the highly relativistic nonequatorial circular motion with 
r = 3M(1 — ^i)~^, cos^ = S2, <p = canst, S\ = 
const ^0, 5*2 = const ^0, 83 = 0. That is, Eq. (43) 
shows the possibility of the spinning particle solution that 
orbit permanently above or below the equatorial plane, 
however, with the small value \6 — 90° | only. One can 
verify that due to the small values Si , 62 solution (43) is 
an approximate solution of the rigorous MPD equations. 

Just as in the case of Schwarzschild's field, the set of 
Eqs. (39) with the corresponding expressions for Ai can 
be integrated numerically in Kerr's field. Fig. 10, as 
an analogy of Fig. 2, shows the plots of r(s) for some 
values of the inclination angle at a = M. Fig. 11 shows 
the dependence of the amplitude and period of the 9- 
oscillation on the Kerr parameter a. The main features 
of Figs. 4-9 are peculiar to the corresponding plots for 
Kerr's field, that are not presented here for brevity. 

According to Figs. 2, 6, 10 when the inclination an- 
gle is equal to 90°, i.e., when spin lies in the equatorial 
plane, so that spin-gravity coupling is equal to 0, the 
corresponding plots coincide with the geodesies. It is an 
evidence of the correct transition from solutions of the 
MPD equations to geodesies. 

Fig. 12 illustrates some typical cases of the equatorial 
motions of a particle with fixed initial values of its coor- 
dinates and velocity but with different absolute value of 




FIG. 2. Radial coordinate vs. proper time for the inclination 

angle 0° (horizontal line r = 3M), 1° (dash and dot line), 10° 
(dash line), and 90° (sohd Unc) at a = 0, £0 = lO"*. The dot 
line corresponds to tiic geodesic motion with the same initial 
values of the coordinates and velocity. 
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FIG. 3. Graphs of the angle 6 vs. proper time for the incli- 
nation angle 0° (horizontal line 9 = 90°), 1° (solid line), 10° 
(dash line), and 90° (dash and dot line) at a = 0, eo = 10~*. 
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FIG. 4. Graphs of the angle vs. proper time at a = 0, 
£o = 10~* for different values of the inclination angle practi- 
cally coincide with the corresponding geodesic plot. The same 
feature takes place for the corresponding graphs t vs. s that 
are not presented here for brevity. 
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0.010- 



0.006- 



0.002- 
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FIG. 7. Radial coordinate vs. proper time for the inclina- 
tion angle 0° (dash and dot line), 1° (dash line), and 10° 
(solid line) at a = 0, eo = 10""', and the nonzero radial ve- 
locity dr/ds w 3.858 x 10~^. The dot line corresponds to the 
geodesic motion with the same initial values of the coordinates 
and velocity. 




s (M) 



FIG. 5. Graphs of Si vs. proper time for the inclination angle 
1° (dash line), 10° (dash and dot line), and 90° (solid line) at 
a = 0, £0 = 10~*. 



FIG. 8. Graphs of the angle 9 vs. proper time for the inclina- 
tion angle 0° (horizontal line), 1° (dash line) and 10° (solid 
line), at a = 0, £0 = 10"*, and dr/ds w 3.858 x 10"''. 




FIG. 6. Trajectories of the spinning particle in the polar 
coordinates for the inclination angle 0° (circle r = 3M), 1° 
(dash line), 10° (dash and dot line), and 90° (solid line) at 
o = 0, £o = 10^*. The dot line corresponds to the geodesic 
motion with the same initial values of the coordinates and 
velocity. The circle r = 2M corresponds to the horizon line. 



FIG. 9. Trajectories of the spinning particle in the polar 

coordinates for the inclination angle 0° (dash and dot line), 
1° (dash line), and 10° (solid line) at a = 0, £o = lO"**, 
and dr/ds ~ 3.858 x lO"'^. The dot line corresponds to the 
geodesic motion with the same initial values of the coordinates 
and velocity. 



9 




3.85- 



3.80-[ ^ ': \ ' 

0.1 0.2 0.3 0.4 0.5 

5 (M) 

FIG. 10. Radial coordinate vs. proper time for the inclination 
angle 0° (horizontal line r = 4M), 1° (dash and dot line), 10° 
(dash line), and 90° (sohd hnc) at a = M, eo = 10"". The 
dot line corresponds to the geodesic motion with the same 
initial values of the coordinates and velocity. 




s (M) 

FIG. 11. Graphs of the angle 9 vs. proper time for the incli- 
nation angle 1° at a = (dash line), a = 0.5M (dash and dot 
line), a = M (solid line); eo = 10""*. 




3.80- ".\ ■ 



0.05 0.10 0.15 0.20 0.25 
s (M) 

FIG. 12. Radial coordinate vs. proper time for the equato- 
rial motions in Kerr's field at a = M with different absolute 
value of spin and the same (common) initial values of the co- 
ordinates and velocity. The circular orbit with r = AM at 

£o ~ 10~" is shown by the long dash line. The solid, dash 
and dot, and dash lines describe the cases when eo is equal to 
0.9 X 10"", 0.6 X 10"*, and 0.2 x 10"* correspondingly. The 
dot line represents the geodesic motion with the same initial 
values of the coordinates and velocity. 



Finally, we remark on a simple conclusion following 
from the nongeodesic curves of a spinning particle pre- 
sented, in part, in Figs. 2-12. Let us consider any point 
on the trajectory of a spinning particle corresponding to 
its proper time Si > (we recall that all curves in Figs. 
2-12 begin at s = 0). Then the geodesic curve can be 
calculated which starts just at this point with the veloc- 
ity that is equal to the velocity of the spinning particle 
at the same point. Also, it is not difficult to estimate the 
deviation of the pointed out nongeodesic curve from this 
geodesic. In principle, it means that in our comparison 
of the corresponding geodesic and nongeodesic curves we 
are not restricted to the trajectories of a spinning par- 
ticle which start in the small neighborhood of the value 
r — r^iJ only (naturally, here we are restricted to the do- 
main of validity of the linear spin approximation) . This 
conclusion may be useful for the generalization of the re- 
sults obtained in this Sec. on other motions of a spinning 
particle. 

V. CONCLUSIONS 

In this paper, using the linear spin approximation of 
the MPD equations, we have studied the significantly 
nongeodesic highly relativistic motions of a spinning par- 
ticles starting near r = r^^^ in Kerr's field. Some of these 
motions, namely circular, are described by the analyti- 
cal relationships following directly from MPD equations 
in the Boyer-Lindquist coordinates. Others, noncircular 
and nonequatorial, are calculated numerically. For real- 
ization of these motions the spinning particle must pos- 
sess the obital velocity corresponding to the relativistic 
Lorentz 7- factor of order l/\/e. All considered cases of 
the spinning particle motion arc within the framework of 
validity of the test-particle approximation when e ^ 1 . 

The situation with a macroscopic test particle moving 
relative to a massive body with 7^ » 1 is not realis- 
tic. However, the highly relativistic values of the Lorentz 
7-factor are usual in astrophysics for the elementary par- 
ticles. For example, if M is equal to three of the Sun's 
mass (as for a black hole), then eq for an electron is of 
order 0.4 x 10"^^ and from Eq. (28) we have the 7-factor 
of order 2 x 10^. Similarly, for a neutrino with the mass 
w ley we have £0 ~ 2 X 10"", 7 « 3 x 10^. 

We can expect the effects of the significant space sep- 
aration of some highly relativistic particles with differ- 
ent orientation of spin. Indeed, the effects considered in 
this paper exhibit the strong repulsive action of the spin- 
gravity interaction. For another correlation of signs of 
the spin and the particle's orbital velocity this interac- 
tions acts as an attractive force. In general, the last case 
is beyond the validity of the linear spin approximation. 
The nonlinear spin effects will be investigated in another 
paper. Also, we plan to show that according to the MPD 
equations significantly nongeodesic orbits of a highly rela- 
tivistic spinning particle, with the 7-factor of order 
exist for the much wider space region of the initial val- 
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ues of the particle's coordinates in a Kerr spacetimc than 
the orbits considered above. However, the corresponding 
calculations are very complicated because this result fol- 
lows from the rigorous MPD equations (1), (2) only, and 



is not common for the approximate equations (2), (7). 

Naturally, it would be interesting to study the possible 
role of the highly relativistic spin-gravity interaction in 
the jet formation. 
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